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1 Introduction 



In the recent study of the Calogero- Sutherland model (CSM)[1, 2] Yangian 
symmetry and Jack symmetric functions have played a central role [3~ 10]. 

Yangians are quantum groups associated with the rational solution of the 
Yang-Baxter equation (YBE)[11]. Realizing the Yangian Y^gl?) in terms of 
the degenerate affine Hecke algebra, Bernard et al. have formulated the spin 
generalization of the CSM[4]. This allows us to understand the integrability 
of the CSM based on the YBE analogously to other integrable models such 
as Heisenberg spin chains and to construct the space of states based on the 
representation of the Yangian[12]. 

The Jack symmetric function is a multi-variable orthogonal symmet- 
ric polynomial. It has been shown that the space of states of the CSM is 
spanned by the Jack symmetric functions[8, 7, 9, 10]. Furthermore, the or- 
thogonality property of the functions allows one to calculate the dynamical 
correlation functions[9, 10]. Remarkably, the resultant density-density cor- 
relation function turns out to provide a certain universal correlation function 
in the random matrix theory at the special coupling constants g = 1/2,1 
and 2[13]. 

The CSM has its one-parameter deformation known as the trigonomet- 
ric Ruijsenaars- Schneider model (RSM)[14]. In the recent paper[15], we 
have considered this model and shown that its spin generalization possesses 
the quantum affine symmetry generated by the level-0 action of U q (gl 2 )[16]. 
The algebra U q (gl 2 ) at level-0 is known as a trigonometric extension of the 
Yangian Y(gl2). Hence one may expect a similar construction of the space 
of states based on the representation of U q (gl 2 ) at the level-0 to the CSM 
done by the Yangian representation. We also have shown that the (spin- 
less) trigonometric RSM is diagonalized by using the Macdonald symmetric 
functions [17], in the same way as the CSM is by the Jack symmetric func- 
tions. The Macdonald symmetric function is known to be a one-parameter 
deformation of the Jack symmetric function. 

The RSM has been discussed by many people in connection to wide 
variety of fields such as sine-Gordon theory[14, 18, 19], Toda theory[14, 
20, 21], G/G gauged Wess-Zumino-Witten (WZW) model[22] as well as 4d 
supersymmetric Yang- Mills theory[23, 24]. 

In this paper, we consider the calculation of dynamical correlation func- 
tions of the trigonometric RSM. In [15], we have briefly reported the results 
on the dynamical density- density correlation function and the one-particle 
retarded Green function in the finite system with particle number N and 
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circumference L, on which particles are confined. We present here the de- 
tails of the calculation. Our method is based on the Macdonald symmetric 
functions and is completely parallel to the one done in the CSM based on the 
Jack symmetric functions[9, 10]. For this purpose, we prepare some formu- 
lae of the Macdonald symmetric functions in section 3. We next investigate 
the thermodynamic limit, N, L — ► oo with N/L fixed, of the dynamical cor- 
relation functions and obtain natural one-parameter deformed expressions 
of those of the CSM. Physical significance of these results is not so clear 
at this stage. However the argument by Gorsky and Nekrasov[22], which 
shows the equivalence between the trigonometric RSM and the G/G gauged 
WZW model with a certain Wilson line insertion suggests that our results 
should be connected with some amplitudes in the latter model. 

We also investigate the statistics of the excitations of the model. We 
show that for a rational coupling constant g, they obey the fractional ex- 
clusion statistics a la Haldane with statistical parameter #[25]. We stress 
the role of the dual description of the states associated with a conjugate 
partition A' and a coupling constant 1/g. Its consistency with a description 
associated with a partition A and the coupling constant g selects the physical 
processes. 

On the other hand, it has been shown by Kawakami and Yang that 
the low energy behavior of the CSM is described by the C = 1 Gaussian 
CFT[26]. However the critical exponents obtained by Ha[10] based on his re- 
sult of the one-particle Green function seem not to agree with those in [26]. 
This is due to the fact that Kawakami and Yang considered the correla- 
tion functions of operators associated with real fermions or bosons, whereas 
Ha treated those associated with the pseudo-particles. Although Ha has 
explained his results based on the low-energy effective theory of the pseudo- 
particles which indicates the fractional exchange statistics, their connection 
to the standard CFT argument is still unknown. Here we settle this point. 
For this purpose, we derive a new selection rule (6.9) for pseudo-particles 
obeying the generalized Pauli exclusion principle, which restrict the possi- 
ble quantum numbers labeling the critical exponents. Under this rule, we 
show that the critical exponents of the Green function both in the CSM 
and in the trigonometric RSM agree with the CFT predictions. We thus 
conclude that the trigonometric RSM belongs to the same universality class 
as the CSM, i.e. the C = 1 Gaussian theory with the compactification ra- 
dius 1/yfg. As discussed in [15], this indicates that the model behaves as a 
Tomonaga-Luttinger liquid[27, 26] in the low energy region. 

This paper is organized as follows. In section 2, we review the trigono- 
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metric RSM and fix the notations. In section 3, we give a brief review of the 
theory of Macdonald symmetric functions and present some useful formulae. 
In section 4, we discuss the diagonalization of the trigonometric RSM and 
investigate the statistics of the excitations for a rational coupling constant 
g. We show that the excitations of the model obey the fractional exclusion 
statistics with statistical parameter g. Section 5 is devoted to a detailed 
description of the calculation of the dynamical correlation functions. The 
thermodynamic limits of the results and their low-energy asymptotic behav- 
ior are derived. In section 6, we consider the finite-size scaling corrections 
to the excitation energy and momentum spectrum and identify the CFT 
describing the low energy behavior of the model. We also derive a new 
selection rule of the quantum numbers. 



2 The trigonometric Ruijsenaars-Schneider model 

The trigonometric Ruijsenaars-Schneider model (TRSM) is an integrable 
quantum mechanical system of N — relativistic particles [14]. Let 0j (j = 
1,2,..,N) be the rapidity variables and Xj be their canonically conjugate 
variables. We impose the canonical commutation relations 

[x j ,0,] = i6j,h (2-1) 

with Ti = 1. We use the representation d\ = —idjdx\. The model is described 
by the following Hamiltonian H and momentum operator P 

H = ^(H_ 1 + H 1 ), P= C -( H _ 1 -H 1 ) (2.2) 

as well as the following N — independent integrals of motion H k (or H-k) (k = 
1,2,..,JV) 

H±k = Yi n h± ( x < - ^) i/2 e- i/cE ^° n h± (*m - xtfixv 

/c{i,2,..,iv} iei iei 

\I\ = k jgl mil 

s sin^fa;,' — x;, T iqlc) 
h ± x t - Xj ) = . ' J+ 2.4 
sm^j;, - Xj) 

Here c being the speed of light, L £ R>o and g £ Q being the coupling 
constant. The particles are confined on the ring of length L in the space of 
{xj}. We normalize the mass equal to one. Notice that H±k (k = 1, 2, N) 
commute with each other 

[H k ,H l ] = U Vk,le{-N,-N + 1,...,N}. (2.5) 
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The Lorentz boost generator B is given by 

1 N 

5 = —$>,■. (2-6) 

c 1=1 

The model possesses Poincare invariance in the sense that the operators 

H, P, and B satisfy the Poincare algebra 1 : 

[H,P] = 0, [H,B] = iP, [P,B] = i^. (2.7) 
In the 'non-relativistic' limit c — ► oo, we recover the Hamiltonian of CSM 
lim(I - Nc 2 ) 

A 1 / d x2 f s „ (tt/L) 2 , , 

j = l l<j<fc<JV L v •? K ^ 

It is also known that the integrals of motions has a deep connection 
with the Macdonald operators. The Macdonald operators Dk(p,t) (k = 

I, 2,..,iV) are defined by[17] 

D k (p,t)=t k w n^^n/ 8 ( 2 - 9 ) 

/c{i,2,..,iv} iei Zt Z i iel 
\I\=k Jgl 

for iV-variables zj (j = 1,2,..,N) with t,p G C x . p® 3 denotes the shift 
operator (p^ J /)(^i, »,zn) = —,P z ji ■■^ z n)- The operators Dk(p,t) (k = 

1,2, ..,iV) can be simultaneously diagonalized by the Macdonald symmetric 
functions (see §3). 

In order to show the connection between and H^, let us set 

p = e- 2v ' Lc , t = p 9 (2.10) 
zj = e l2 ~r x \ = e ^ Zj ^h. (2.11) 



Then, by using the function 



A = TT ( z 3 1 Zk ' l P) QO (2 i2) 

N ^t, (tZj/Zklp)^ 



J One should note that the Poincare invariance does not guarantee that the Hamiltonian 
(2.2) defines an interacting system of relativistic particles satisfying the causality. In this 
paper, we use the terminology 'relativistic' only in the kinematical sense. 
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(x;p) 00 =l[(l-xp n ), (2.13) 

n=0 



tT k(N-i)/2 Dk ( p ±i^±v- ) k=l,2,..,N. (2.14) 

3 Macdonald symmetric functions 

In this section we briefly review the Macdonald symmetric functions [17]. 
Let A = (Ai, • • • , Ajv), Ai > ••• > Ajv > 0, Xj 6 Z be a partition. The 
Macdonald symmetric function P\(z;p,t) = P\(zi, ..,zw,p,t) is an eigen 
function of the Macdonald operators, 

D k (p ±1 ,t ±1 )P x (z;p,t) = ( l[t N - l 'p X ')Px(z;p,t). (3.1) 

h< — <ik '=1 

The Macdonald symmetric functions form an orthogonal basis of a ring 
of symmetric polynomials Ajv in the N variables z = (z\, Z2, ..,zn) with 
coefficients in Q(p, /). Let /(A) be the length of A, i.e. /(A) = max{j \Xj > 0}. 
To each partition A, we assign a Young diagram V(X) = < i < 

K^)' 1 — J — hJ ^ Z>o}. We denote the conjugate partition of A by 

A' corresponding to the transpose of the Young diagram T>(X). For each cell 
-y = of V(X), we define the quantities (2(7) = A; — j, 0/(7) = j — 1, 

/( 7 ) = A$ - i and Z'( 7 ) = i - 1. 

Define the scaler product in Ajv by 

1 f N d 

<f,g>' P , tiN = ^Jn^m^N 

= —(constant term in f(z)g(z)A]\f), (3-2) 



with f(z) = f(l/zi, .., 1/zn), then the orthogonality of the two Macdonald 
symmetric functions is stated as follows. 

< P„ P A >; jt = A N N(X), (3.3) 



where 



one has [28] 

A -1 / 2 77 A 1 / 2 - 
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where 

A N = < 1, 1 >' P ,l;N = ft (3.4) 

/. A (p,o= n^-^ 7 ^ 1 )' ( 3 - 5 ) 

„ 1 _ v *'(l)t N - l '(l) 

^) = n 1 _ p X) + i^ ) -i - ^ 

The P\(z;p,t) satisfies the following duality property: 

^2P x (z; P ,t)P xl (w;t,p) = l[(l + ZiW j ). (3.8) 

In the calculation of dynamical correlation functions (§5), the following 
expansion formulae[17, 29] are essential. 

f>f = (i-P m ) £ ^M p a(.;p,0, (3-9) 



|A| = 



where a G C x and 



A i(«^;p)oo ^ h X '(t,p) 



x\ P ,t)= n (i'' w -p a ' w ), (3.H) 

76A 
7^(1,1) 

(4^ = n(^' (7) -P a ' (7) «)- (3-12) 

7GA 

4 Diagonalization 

We here consider the diagonalization of the TRSM and investigate the statis- 
tics of the elementary excitations. 

From (2.14) and (3.1), one can diagonalize H and P by the functions 

$jv ;A (z) = A]fp A (z;p,i) (4.1) 
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with the following eigenvalues 



N n N 



E N (\) = c 2 ]Tch-f, P N (\) = c]Tshf , (4.2) 

where 



. c " c 



«i = T^ + »(^T 2 -•*)}• < 4 ' 3 » 

The model thus can be regarded as an ideal gas of N — relativistic pseudo- 
particles with the pseudo- rapidities (4.3). 

The formula (4.3) implies the minimum difference between two consec- 
utive pseudo-rapidities 

Oi - 6 i+ i > ^g- (4.4) 

When g = there are no restrictions, and more than two pseudo-particles 
can occupy the same quantum sate so that the statistics is bosonic. In the 
case g = 1, a pseudo-particle occupy one of the quantized states by the unit 
2tt/L obeying Pauli exclusion principle so that the statistics is fermionic. 
For generic g, the pseudo-particles can be said to obey the generalized Pauli 
exclusion principle[25]. 

1/2 

The ground state is given by the function ^>q{z) = corresponding 
to the empty partition A = 0. Note Pg)(z;p,t) = 1. The ground state 
energy and momentum eigenvalues are obtained as 

tp(o) 2 ^gN irg 

L N ' = c sh— — /sh— (4.5) 
Lc Lc 

and pffi = 0, respectively. One should note that the ground state can be 
described as a filled Fermi sea with pseudo-momenta P- ^ = c sh^f with 
— Op < 0J <0 F j = l,2,..,iV, where 

e F = 7rg(iV - 1) . (4.6) 

Let us next consider the excitations. We consider the case with the 
rational coupling g = r/s with r and s (j^ 0) being coprime integers. 

In order to describe the excitations, we assign a motif to each set of 
pseudo-rapidities {9j} j = 1,2,..,N. A motif is a sequence of numbers or 
1 with totally N of Is. It can be constructed as foilows[10]. Consider a one- 
dimensional lattice whose lattice spacing is 2ir/sL, and assign each lattice 



8 



point with 1 if that point coincides with the value of one of the pseudo- 
rapidities, otherwise assign 0. The allowed number of 0s between each pair 
of Is is r — 1 + ns with n G Z>o- Among the r — 1 + ns 0s, we regard r — 1 
consecutive 0s as to be bounded to each 1 and the other 0s as unbounded. 

We define a quasi-particle as 1 which lattice point 9 satisfies 9 > 9p or 
9 < — 9p. Out of r — 1 + ns 0s, we define a quasi-hole as s consecutive 
un-bounded 0s whose lattice points are lying in the region [—9p,9p]. 

By using the duality property (3.8), one has a dual description of the 
excited states, where the roles of quasi-particles and -holes are exchanged. 
More precisely, let p = t,t = p. Then / = p 1 ^ 9 and D k (p,t) = D k (t,p). 
The operator Dk(p,t) are diagonalized by the dual Macdonald polynomials 
P\i(t,p) yielding the eigenvalues Y^f=i Ylt=i p N ~ tl t X ' . Hence if one defines 
the dual integrals of motion H±k by 

H ±k = t^ N -^Af Dk (pJ)A-^ (4.7) 

with Ajv = Ajvlp^t, then for H = y(_ff_i + Hi) one has 

N cj 

HAfPyipJ) = (/5 2 ^ch^)A^ /2 P v (p,0 (4.8) 

3 = 1 C 

with ~c = c/g and ^ = ^{\>. + - j))- 

For the dual rapidities {f^j}, we define a dual motif in the same way as 
in the above except for changing r and s. In the dual picture, the quasi-hole 
is defined as 1 which lattice point $ satisfies $ > -dp or $ < — -dp and the 
quasi-particle is defined as r consecutive un-bounded 0s lying in [—-dp, -dp] 
with 'dp = 

We require that a physical state should be described by the above two 
descriptions consistently. Namely, for a given A, we regard a correspond- 
ing state as physical only if the number of allowed quasi-particles and -holes 
obtained in the one description coincides with those obtained in the dual de- 
scription. For example, in the case g > 1, the partition A = (r — 1, .., r — 1) 

s 

seems to be able to give excitations of r — 1 quasi-particles and s quasi- 
holes. However in the dual picture, this case is forbidden because of the 
condition that the un-bounded 0s corresponding to quasi-particles should 
be in [—-dp, -dp]. In such case, we regard the excitations are forbidden. 

Now require that the number of quasi-holes should be an integer. Then 
from the above requirement and (4.3) one can find that for quasi-holes exci- 
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tation, there is the minimum number, which is equal to r, and it is accom- 
panied by s quasi-particles excitation. The minimal partition corresponding 
to this excitation is given by (r, .., r). This indicates that the quasi-particles 

s 

and -holes obey the fractional exclusion statistics with the statistical pa- 
rameter g = r/s as in the CSM[25, 10]. We will show later (§5) that the 
resultant physical state under this requirement is consistent to the descrip- 
tion of the intermediate state both in the density-density correlation function 
and one-particle Green function. 

If one allows the change of the total particle number N, there is an- 
other way of obtaining excitations. Namely, quasi-holes can be excited by 
destroying some of the pseudo-particles in the ground state configuration. 
As discussed by Ha in the CSM[10], this case also leads to the same frac- 
tional exclusion statistics as in the above. Namely, for integer number of 
quasi-holes excitation, the quasi-hole should be excited by the unit con- 
sisting of r quasi-holes. The one unit excitation is then possible by the s 
pseudo-particles destruction from the ground state. 

In general, allowing the change of the total particle number, excitations 
are obtained by mixing the above two types of excitations. A typical exam- 
ple is the intermediate state of the one-particle Green function, where one 
pseudo-particle destruction from the ground state configuration should be 
taken into account (see §5). 

One should also note that the pseudo-rapidities (4.3) obey the following 
Bethe ansatz like equations: 

N 

L9j = 2-Klj + Tr(g sgn(0j - 0i), (4.9) 

;=i 

with Ij = Xj + —j. In the thermodynamic limit, (4.9) yields the relation 

1 = g PP (6) + PH (6), (4.10) 

where pp{9) (resp. />#(#)) are the particle (resp. hole) density. This also 
indicates the fractional statistics, i.e. one unit increase of the density of 
the particles of rapidity 9 is accompanied by decrease of the hole density by 
g— unites[30]. 

5 Dynamical correlation functions 

The orthogonality property of the Macdonald symmetric functions allows us 

to calculate the dynamical density-density correlation function (0|p(£, 7")p(0, 0)|0) 
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as well as the one-particle retarded Green function (0|$t(£, r)$(0, 0)|0), 
where £ and r are real space and time coordinates, respectively. 



5.1 Finite system 

The dynamical density-density correlation function is defined by 

«*o.o)|o) = *<Mp* (5.D 

JV\U|U)jv 

Here the density operator p(£,r) is defined by 

p(£,T) = e^- p «V(0,0)e-^- p «, 

N 

p(0,0) = ^6( Xj )-N/L, (5.2) 

1/2 

and |0)j\r denotes the iV— particle ground state, i.e. |0)jv = . Combining 
(3.9) with the formula 

-, oo JV 

p(o,o) = I EE(^ m + ^ ( 5 - 3 ) 

m = l j=l 

one can expand the state p(0,0)|0)jv in terms of the N— particle eigen states 
&N;\( Z ) (4-1)- Then using the orthogonality (3.3), we obtain the density- 
density correlation function as follows[15]. 

(o|Ke,r)Ko,o)|o) 

_ 2 (l-p|A|)2 (x A (jM))2 £n{x)t 



(5.4) 

where £jv(A) = En(X) — E$ and Pjv(A) are the excitation energy and 
momentum, respectively. 

One should remark that the factor % A (p,i) in (5.4) vanishes if the Young 
diagram V(X) contains the lattice point (s + l,r + 1). Because of the argu- 
ments in §4, this indicates that only the excited states which can contribute 
to the intermediate state are those consist of minimal s — quasi-particle and 
r — quasi-hole excitations. As in the CSM, this supports the fact that the 
quasi-particles and the quasi-holes obey the fractional exclusion statistics 
with statistical parameter g = r/s. 
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The one-particle retarded Green function is defined by 




Here $t and $ are the creation and annihilation operators of the pseudo- 
particles, respectively, and 



tft (fjT ) = e ^-^%t( ,0) e -^- p « (5.6) 



The action of $(0,0) on the ground state |0)jv+i is defined by 



*(0,0)|0)jv+i = A 




(1/^;p)oo(^;p)oq y/ 2 A i/2 



Z N+1 =1 



(5.7) 



and 



JV +i<0|* t (0,0) = (*(0,0)|0)jv+i). 



(5.8) 



By using (3.10), one can expand the state (5.7) by $n-,\( z ) again- The result 
is 



The factor (i _1 )|f' in (5.10) vanishes if the Young diagram V(X) contains 
the lattice point (s, r + 1). This indicates that only the states which contains 
minimal s— 1 — quasi-particle and r — quasi-hole excitations can contribute to 
the intermediate state. Since $(0,0) annihilates one pseudo-particle from 
the N + 1— particle ground state, this feature is again consistent with the 
fractional exclusion statistics with statistical paprmeter g = r/s discussed 
in §4. 



tf(0,0)|0)jv +1 




(5.9) 



Then by (3.3), we obtain[15] 

(0|Wt(e,r)*(0,0)|0> 




(5.10) 
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5.2 Thermodynamic limit 

Let us next consider the thermodynamic limit N,L — ► oo, with N/L = n 
fixed. Let us first consider the dynamical density-density correlation func- 
tion. Noting the remark below (5.4) and using the formulae given in §A.l, 
we obtain the following expression. 

(0|Ke,r)K0,0)|0> = C(n/ da t )(H dfr) sh 2 — (£ a t + £ fr) 

2=1 J = l 2=1 J = l 

xF{s,r;g\{ai,^})co5{Vr,s0e- i£T '- T . (5.11) 

Here we set a 4 - = A 4 -/iV and (3j = Aj/iV. The quantities 7-V )S and £ rjS are the 
thermodynamic limit of Pjv(A) and £jv(A), respectively: 

;=i j=i 

^ = c 2 E^w + c 2 E e ^0') ( 5 - 13 ) 
i=i j=i 

with 

p P (i) = 2sh— -ch — v g » ( 5 - 14 ) 

2 , irnqfl; , irnq(l — 8;.) 

PhU) = -sh -^ch — ^ ^, 5.15 

g c c 

ep (Q = 2sh^sh 7rra(a ; + g) , (5.16) 

^(j) = -sh^sh 7rrag(1 -^ ) . (5.17) 
5 c c 

Due to the argument below (5.4), pp and ep (resp.p# and e#) are regarded 
as the quasi-particle (resp.quasi-hole) contribution to the excitation momen- 
tum and energy. The normalization constant C is obtained as 

C=4(^)%^-W(rM( S ,r; 5 ), (5.18) 

, x r'(o)L fc (i/fl) 

A(l,k;g) = —, ^ — K T 1JLl . (5.19) 
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The form factor F(l,k,g\{pi,Tj}) is evaluated as 
F(l,k-g\{a % J 3 }) 

= riri sh " 2 -(««+^) 

3=1 i=l C 

nU,- |sh ^(a,- - a 3 T Iti<i l sh ^ ~ W 9 
n! =1 [sh ^a,-sh ^(a t - + g)Y-9 ^[sh ^ ~ ft)] 1 " 1 ' 3 ' 

(5.20) 

Similarly, by using the formulae in §A.2, the thermodynamic limit of the 
one-particle Green function is obtained as 



(01^,^(0,0)10) 



s-1 



=v(n r^)(n / d^eM { ^ (e»«+e^)> 

2=1 J = l 2=1 J=l 

xF(s - l,r-g\{a t J 3 })e-^-^- v ^\ (5.21) 

where 

I i\ i / 7TTl \ r + S — 1 2jrno -. Jrng(g — 1) 

X> = 4/( s - 1 )+^(— ) (1 - e -^-i e ^^-^_U( s _ fl X5.22) 

It should be remarked that in the non-relativistic limit c — ► oo, the 
results (5.11) and (5.21) coincide with those obtained by Ha in the CSM[10] 
up to scalar multiple. 

5.3 Low energy asymptotic behavior 

The exact results obtained in §5.2 allows us to analyze their low energy 
asymptotic behavior. Here 'low energy asymptotic behavior' means that 
we analytically continue the time variable r to the imaginary time —it 
and consider the leading behavior for r 1. The leading contribution 
is obtained from the saddle points, i.e. the points satisfying 8 T)S = (resp. 
S r<s -i = 0) for the density- density correlation function (resp. for one-particle 
Green function). From (5.13), there are r + 1 such points, where all of a 4 - and 
r — to (to = 0, l,..,r) of (3j are and remaining to of (3j are 1. Expanding 
around these points, we obtain 

/ k — m 



V k ,ii ± iS k ,ir = 2-irn a t + ^ /3j] 



i=i j=i 
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-27rn(± ^ Lo a + 2mpF^ + higer terms, (5.23) 



a = l 

where we set (± = £ch ^ ± icr sh w a = 1 - (i a and p F = |sh ^ 
denotes the fermi momentum. The expansion of the dynamical correlation 
functions (5.11) and (5.21) can be carried out in the same way as in the 
CSM[10]. We obtain their leading order behaviors as follows. 

(o|Ke,r)Ko,o)|o) 

« C If (l|,,r){(^) 2 + (^) 2 } + £ C m (^_) m2/fl cos(p F mO, 

(5.24) 

(o\&(Z,T)y(o,o)\o) 

\— / 1 \ {m— g) 2 /g ( 1 \' m2 /g 
~ Z^M^^C+J 1 2^1 J expi(p F m£ + i/zr). (5.25) 

Here fi = c 2 ch being the chemical potential and 

c = 9 .<-» c( =)-(^) 2 (M6 ) 



C m =,— j '^/sh^j '-I?(l| a ,r-m)If(l|m), 

(5.27) 



Xlf(l|s- l,r -m)If(l|m) (5.28) 

where 

if & - m) = (n r ( n r n ff (<*,■ + 

i=l 1/0 j=l 1/0 i=l 3 = 1 



nU«!- fl nfcr/3j" 1/fl 

/ k — ra 

X exp{-z(^ a, + ^ /3j)} X X m (a,-, /%), (5.29) 
i=i j=i 

w= (n r n tl Uh ~^i' 9 e ^- z £ ^ 

a^ lT=l^a a = l 
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In (5.29), the function I m (a 8 ,/3j 



( sh 2 ^ for m ^ 

The factor lf(l|s — l,r — m) in (5.28) is obtained from lf(l\l,k — m) by 
setting / = s — 1, k = r and replacing I m (a 8 ,/3j) with 

| expj^l-^XErV' + £^0} form=0; 
j exp|^p(l- 5 )} form t^O 

In the derivation of (5.24) and (5.25), we carried out analytic continuation [31, 
10] 



/ 1 \ Ql 2 — 2lm + 2S m r +m 2 la „ 

X((z\l,m) = [-J If(l|/,m), (5.31) 

l 3 2 (z\m) = (-) m2/S lf(l|m). (5.32) 



The asymptotic behaviors (5.24) and (5.25) coincides with those of the 
CSM[10]. 



6 Finite-size scaling 

The exact energy spectrum obtained in §4 allows one to analyze the finite- 
size scaling of the model in the thermodynamic limit N, L — ► oo and N/L = 
n being fixed[32, 33, 34, 26, 35]. From this we identify the CFT describing 
the low-energy behavior of the model. 

First of all, from (4.5) we obtain the finite-size correction to the ground 
state energy as 

lim4 0) = J L eo -^ + 0(-l), (6.1) 

3 

where we identified the ground state energy density with £o = — sh and 

Trg C 

the velocity of the elementary excitation with v = c sh 

Comparing the result (6.1) with the general theory of the finite-size 
scaling[33, 34], one may suspect that the central charge is given by g. How- 
ever this is not the correct identification[26]. The central charge should be 
identified with 1. This can be justified, for example, by calculating the low 
temperature expansion of the free energy from (4.9). Instead of doing this, 
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we here justify it by deriving the whole conformal dimensions associated 
with the elementary excitations. 

Let us consider the elementary excitations associated with the particle 
number change by AN and the transfer of the AD particles from the left 
Fermi point to the right one. We evaluate the following finite-size correction 
to the excitation energy and momentum. 

r N +* N 2vr/ / N + AN + l\ \ (0 )1 

l im [ c 2 J2 ch — (g(j- ^ ^ )+AD)-E^\ 



2lrV 19 A l\l2 i 1 A „2l , 1 



= f iAN+— [±AN' + -AD'\ +0(jz ), (6.2) 

2vrch ^ 1 
= 2p F A£ + ^^AiVA£ + 0(— ), (6.3) 

where fi = c 2 ch is the chemical potential. Adding the energy 2j ^ L (N + + 

27rch 11311 _|_ 

iV~) and the momentum L c (N + — N~ ) of the quasi-particles and quasi- 

holes exciting near the Fermi surface, we obtain the finite-size corrections 
up to the order 1/L: 

AE = nAN + ^[^AN 2 + ^AD 2 + N + + N~}, (6.4) 

27Tctl 7T ^ n 

AP = 2tt Pf AD + j^^ANAD + N + - N~^. (6.5) 

From these, we obtain the right and left conformal dimensions as[15] 

h ± (AN; AD; JV±) = ^ [ff ± \n ± , (6.6) 

where we set R = 1/y/g. The conformal dimensions with = are 
nothing but those of the U(l) — primary fields in the C = 1 Gaussian theory. 
We hence identify the desired CFT with the C = 1 Gaussian theory with 
compactification radius 1/^/g. 

The standard theory of CFT then yields that correlation functions of 
any operator O can be expanded by those of the primary fields and the 
descendent fields as follows. 

e t PF AD( 

< 0(£,O0(O,O) >« £ const. h+ h . (6.7) 
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It is instructive to compare this with the exact results obtained in §5. In 
the case 0(£,t) = p(£,t), we have no particle number changes so AN = 0. 
Then the formula (6.7) with (6.6) and the identification AD = m agrees 
with the low-energy behavior of the dynamical density-density correlation 
function (5.24). 

On the other hand in the case of the one-particle Green function, we 
have particle number change AN = 1. In such a case, we has to consider 
the statistics of the pseudo-particle annihilated by Such statistics is 
formulated as the boundary condition of the wave function[26, 35]. We 
assume that the wave function of the pseudo-particles with the rapidity 
(4.3) obey the periodic boundary condition 

e l6 i L = 1. (6.8) 

This yields the following selection rule 

AD = |AJV (mod 1). (6.9) 

Substituting (6.9) with AN = 1 into (6.6), we obtain 

(AD + g) 2 , AD 2 

This agrees with the exact result (5.25), if one sets AD = m. We hence 
conclude that the low-energy effective theory of the TRSM (as well as CSM) 
is described by the C = 1 Gaussian theory with compactification radius 
1/y/g and selection rule (6.9). 

The meaning of the selection rule (6.9) in the context of the Gaussian 
theory is as follows. According to the argument in [7], a pseudo-particle 
carries a flux irg. The flux twists the boundary condition of the boson 
field (j)(cr,T) describing the Gaussian theory from <f)(a + 2iri,T) = (f>(cr, r) + 
27ri2m, m G Z to cf)(a + 27ri,r) = <f>(a, r) + 2itR{m + 9 ^ N ). This twisted 
Gaussian theory is similar to the low-energy effective theory discussed by 
Ha[10]. 

7 Conclusion and discussions 

We have discussed a description of the excitations and a calculation of the 
dynamical correlation functions in the trigonometric Ruijsenaars- Schneider 
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model. In the process, we have applied the theory of the Macdonald sym- 
metric functions. We have shown that for the rational coupling constant 
g, the excitations of the model obey the fractional exclusion statistics with 
statistical parameter g. We also have obtained the one-parameter defor- 
mations of the dynamical correlation functions of the Calogero-Sutherland 
model. From the finite-size scaling analysis, we have also shown that the 
trigonometric Ruijsenaars-Schneider model is effectively described by the 
C = 1 Gaussian theory with R = 1/y/g in the low-energy region under a 
new selection rule for the quantum numbers labeling the critical exponents. 

Recently, Gorsky and Nekrasov have shown the equivalence between 
the trigonometric Ruijsenaars-Schneider model and the G/G gauged WZW 
model with a Wilson line insertion[22]. This is a natural deformation of 
the equivalence between the Calogero-Sutherland model and the 2d Yang- 
Mills theory[36, 37]. To show this equivalence at the amplitude level should 
deepen the understanding of both the G/G gauged WZW model and the 
Ruijsenaars-Schneider model. 

In addition, in the classical level, the hyperbolic Ruijsenaars-Schneider 
model is known to have a deep connection with the sine-Gordon theory[18, 
19]. This connection should be examined at the quantum level. 

Another possible problem is a connection to the random systems such 
as random matrix model, 2d disordered system in condensed matter physics 
and the quantum chaos. The latter systems indicate certain universal behav- 
ior and the Calogero-Sutherland model provide these systems with universal 
functions. The parallel structure of the Ruijsenaars-Schneider models to the 
Calogero-Sutherland seems to suggest the existence of some one-parameter 
deformation of the random systems. 

We hope to return to these problems in future. 
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Figure 1: A Young diagram is divided into three sub- diagrams. This diagram 
describes the case g = 3/4. 



memorial foundation. 



A Formulae for the thermodynamic limit 

A.l Dynamical density-density correlation function 

Due to the remark below (5.4), only the partitions which Young diagram 
does not contain the cell (s + 1, r + 1) can contribute to the density-density 
correlation function. Following Ha[10], we divide such Young diagrams into 
the following three sub-diagrams (Fig.l). 

A={(i,j)\l<i<s,l<j<r}, 

B = {(i,j)\l<j<r,s + l<i<X' j }, (A.l) 
C = {(i,j)\l<i<s,r + l<j <Xi}. 

Now let us first consider the factor A/"(A) appearing in (5.4). According 
to the decomposition (A.l), A/"(A) is decomposed into the following three 
factors. 

= nn /.In-, , (a.2) 

j=l t=l F 
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r 3 1 - rJ-ifN-i+l 

= n n / J tN -, 

j = li=s + l L ~ Pt 

r Tt(N _ s + 1 + (j _ i) /g) r f (iv - a;. + j ig) 

f = \ UN - A;- + 1 + (j - l)/g)T t (N - s + j I >g) ' 1 ' ' 

s A 8 1 _ ,_i f N-i+l 

= P ''_;,-, 

1 = 1 j = r-\-l ± 

* T v (\ +g(N-i + l))T p (r + 1 + g(N - i)) 
i\ T P (r + g(N-z + l))T p (\i + l + g (N-i)y [ ' ) 



where 



In the thermodynamic limit N, L — ► oo with iV/i = n fixed, we obtain 

■Af(^) -+ 1. (A.6) 

Noting p e c ^ t e c and the formula 

T t (z + a) fl-t^a 



T t (z) Vl- f 

for |z| — ► oo, we obtained the thermodynamic limit of the remaining factors 
as follows. 

2nng(l-l3 J ) 

m)~fi(^^P) l "-\ (A.8, 

3 = 1 1 - e c 

s _ 2irn(g + ct 8 ) 

^) - n( 1 " e " *l r 1 - ( a - 9 ) 

i=i 1 - e 

Similarly, the factor X A (P? i s given by the product of the following three 
factors 

**m = nnc*''" 1 -^'" 1 ) 

3 = 1 i=l 

(p;p)oo 



^_ys p rs(r-l)/2^ _ t y(s-l) 



(p r+1 ;p)c 
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-A' 

j 



(A.11) 



r K 

x b (p^) = n n (^-p^ 1 ) 

i=l i=s+l 

= (_)"+E j ^ p -"('-i)/2+E; =1 0-i)^ JJ(i _ t y 

r f (A;. - (j - i)/ g ) 
x r t (, - (j - 1)/ 5 ) ' 

x c (p,t) = fl II (< , '- 1 -p , ' _1 ) 

z = l j = r-\-l 
= i -rs( S -l)/2 + X;* =1 («'-l)Ai 

Next we consider the factor h\(p,t) = h^(p,t)hs(p,t)hc(p,t). The factor 
hj[(p,t) is evaluated as 

hA(p.t) = n n(i - p a '-vh +i ). (A.i3) 
j=i i=i 

In order to evaluate h&(p,t), we decompose the diagram £> further into the r- 
cells such that the /— th cell is given by < j < /, AJ +1 + 1 < i < AJ}. 

If AJ = AJ +1 , then the /— th cell is empty. We thus obtain the following 
expression. 

r I K 



h B (p,t) = nn n a-/*-^ 



/=lj=li=A| +1 +l 

f[(l- tyW-^TtW + l-j/g) 
[ fr UK -Xj + l + U-i- 1)1 9) 

x ll UK - a;. + 1 + U - i)/g) " (A " 14) 
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Here note X' r+1 = s. Similarly, we obtain 

S I A; 



h c (p,t) = nn n (i-/ < - j ^- i+i ) (a.is) 

;=ii=i j=A i+1 +i 



n^-f) 



1=1 



8 (A 8 -r) r p( A » -9(i- 1) 

r P (ff) 



.A'. -is 



-A- T p (Xi — Xj + g(j — i)) 
X j-}.T p (X i -X j+ gU-i+l)y (A - 16) 

In the same way, we obtain 

hMt, P ) = f[f[(l - P X '- 1+1 t X '^) (A.17) 

3 = 1 i=l 

hB>(t, P ) = nn ri n-p x '- i+i t } 

l=lj=li=X' l+1 +l 

r Tt(K - A ;. + (J - i)/g) 

x llr t (A:.-A;. + (j - i + i)/,)' (A - 18) 

m*,^ = nn ri (i- P x '- 3+i t x ^) 

1=1 t=l j=A i+1 +l 

= f[a-py {x '- x ' +i %(\ + i-9i) 

i=l 

* T p (X l -X J + l +g (j- i -l)) 
££ T p (X l -X J + l + 9 (j- i )) ■ {AA9) 

Combining (A. 10) ~ (A. 19) and taking the thermodynamic limit, we 
obtain the following results 

h A (p,t)h A ,(t,p) ^ UJ rg i(r) l\l = \ {1 } 



fr T%s - (J - I)/ g) 

x l\Tm-(j-iy 9 y (A - 20) 
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(x B (p,t)f f 2irg^-r{s-i) ^ T(l/g) 



h B (p,t)h B ,(t,p) \LcJ j.J. iT 2(s-(j-l)/g) 

J_ _ 4jrn(j _ 1)/3 27rng/3 

J] e s i (l - e 



x ll e 

r 



hc{p,t)h c <{t,p) \Ld l\ T*(r - </(i - 1)) 

-4:irn(i — l)ga^ 2irna. 

e " 

i=l 



X n e _ e r^)5-i 



x^l-e c fa. (A.22) 

One should note that all the factors 1/L are absorbed into n in (5.18) with 
N appearing from the replacement J2\, ~^ N J da{ and J2\' ~^ N J d(3j. 

A. 2 One-particle Green function 

According to the same decomposition of the Young diagram, we obtain in 
the thermodynamic limit 



hA(p,t)hA'(t,p) K Lc i=1 j =1 



A r 2 (, - (j - l)/ g ) 

x M r=»(i - c/ - ( 3) 



(r 1 )^) 2 . ^2^ r r ( .-i)^ r(i/ 5 ) 



h B (p,t)h B ,(t,p) \LcJ l = l r 2( s - (j - i)/ fl ) 

-4jrn(j-l-g)/3 J 2ixngP J ^ ^ 

j = l 



X n e 3 i (l - e r^) 1 ^- 



x]I( 1 - e ' ( A - 24 ) 

i<3 



((t- 1 )^ ? / 27T X 2r- S (r-l)-l- a T(g) 

h c (p,t)h c ,(t,p) ~* \Lc) fj^ir-gi) 
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s-1 

____ 271-7101, _. 

x ri( i - e — 

i=l 

xl[(l-e ^f 9 . (A.25) 
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